Improved super-resolution ribosome profiling reveals prevalent translation of upstream ORFs and small ORFs in Arabidopsis.

A crucial step in functional genomics is identifying actively translated ORFs and linking them to biological functions. The challenge lies in identifying short ORFs, as their identification is greatly influenced by data quality and depth. Here, we improved the coverage of super-resolution Ribo-seq in Arabidopsis (Arabidopsis thaliana), revealing uncharacterized translation events for nuclear, chloroplastic, and mitochondrial genes. Assisted by a transcriptome assembly, we identified 7,751 unconventional translation events, comprising 6,996 upstream ORFs (uORFs) and 209 downstream ORFs on annotated protein-coding genes, as well as 546 ORFs in presumed noncoding RNAs. Proteomic data confirmed the production of stable proteins from some of these unannotated translation events. We present evidence of active translation from primary transcripts of trans-acting small interfering RNAs (TAS1-4) and microRNAs (pri-MIR163 and pri-MIR169) and periodic ribosome stalling supporting cotranslational decay. Additionally, we developed a method for identifying extremely short uORFs, including 370 minimum uORFs (AUG-stop), and 2,921 tiny uORFs (2 to 10 amino acids) and 681 uORFs that overlap with each other. Remarkably, these short uORFs exhibit strong translational repression as do longer uORFs. We also systematically discovered 594 uORFs regulated by alternative splicing, suggesting widespread isoform-specific translational control. Finally, these prevalent uORFs are associated with numerous important pathways. In summary, our improved Arabidopsis translational landscape provides valuable resources to study gene expression regulation.

Phosphorylation Dynamics in a flg22-Induced, G Protein-Dependent Network Reveals the AtRGS1 Phosphatase.

The microbe-associated molecular pattern flg22 is recognized in a flagellin-sensitive 2-dependent manner in root tip cells. Here, we show a rapid and massive change in protein abundance and phosphorylation state of the Arabidopsis root cell proteome in WT and a mutant deficient in heterotrimeric G-protein-coupled signaling. flg22-induced changes fall on proteins comprising a subset of this proteome, the heterotrimeric G protein interactome, and on highly-populated hubs of the immunity network. Approximately 95% of the phosphorylation changes in the heterotrimeric G-protein interactome depend, at least partially, on a functional G protein complex. One member of this interactome is ATBα, a substrate-recognition subunit of a protein phosphatase 2A complex and an interactor to Arabidopsis thaliana Regulator of G Signaling 1 protein (AtRGS1), a flg22-phosphorylated, 7-transmembrane spanning modulator of the nucleotide-binding state of the core G-protein complex. A null mutation of ATBα strongly increases basal endocytosis of AtRGS1. AtRGS1 steady-state protein level is lower in the atbα mutant in a proteasome-dependent manner. We propose that phosphorylation-dependent endocytosis of AtRGS1 is part of the mechanism to degrade AtRGS1, thus sustaining activation of the heterotrimeric G protein complex required for the regulation of system dynamics in innate immunity. The PP2A(ATBα) complex is a critical regulator of this signaling pathway.

A Practical Guide to Inferring Multi-Omics Networks in Plant Systems.

The inference of gene regulatory networks can reveal molecular connections underlying biological processes and improve our understanding of complex biological phenomena in plants. Many previous network studies have inferred networks using only one type of omics data, such as transcriptomics. However, given more recent work applying multi-omics integration in plant biology, such as combining (phospho)proteomics with transcriptomics, it may be advantageous to integrate multiple omics data types into a comprehensive network prediction. Here, we describe a state-of-the-art approach for integrating multi-omics data with gene regulatory network inference to describe signaling pathways and uncover novel regulators. We detail how to download and process transcriptomics and (phospho)proteomics data for network inference, using an example dataset from the plant hormone signaling field. We provide a step-by-step protocol for inference, visualization, and analysis of an integrative multi-omics network using currently available methods. This chapter serves as an accessible guide for novice and intermediate bioinformaticians to analyze their own datasets and reanalyze published work.

Transcription Factor Dynamics in Cross-Regulation of Plant Hormone Signaling Pathways.

Cross-regulation between hormone signaling pathways is indispensable for plant growth and development. However, the molecular mechanisms by which multiple hormones interact and co-ordinate activity need to be understood. Here, we generated a cross-regulation network explaining how hormone signals are integrated from multiple pathways in etiolated Arabidopsis (Arabidopsis thaliana) seedlings. To do so we comprehensively characterized transcription factor activity during plant hormone responses and reconstructed dynamic transcriptional regulatory models for six hormones; abscisic acid, brassinosteroid, ethylene, jasmonic acid, salicylic acid and strigolactone/karrikin. These models incorporated target data for hundreds of transcription factors and thousands of protein-protein interactions. Each hormone recruited different combinations of transcription factors, a subset of which were shared between hormones. Hub target genes existed within hormone transcriptional networks, exhibiting transcription factor activity themselves. In addition, a group of MITOGEN-ACTIVATED PROTEIN KINASES (MPKs) were identified as potential key points of cross-regulation between multiple hormones. Accordingly, the loss of function of one of these (MPK6) disrupted the global proteome, phosphoproteome and transcriptome during hormone responses. Lastly, we determined that all hormones drive substantial alternative splicing that has distinct effects on the transcriptome compared with differential gene expression, acting in early hormone responses. These results provide a comprehensive understanding of the common features of plant transcriptional regulatory pathways and how cross-regulation between hormones acts upon gene expression.

Brassinosteroids modulate autophagy through phosphorylation of RAPTOR1B by the GSK3-like kinase BIN2 in Arabidopsis.

Macroautophagy/autophagy is a conserved recycling process that maintains cellular homeostasis during environmental stress. Autophagy is negatively regulated by TOR (target of rapamycin), a nutrient-regulated protein kinase that in plants is activated by several phytohormones, leading to increased growth. However, the detailed molecular mechanisms by which TOR integrates autophagy and hormone signaling are poorly understood. Here, we show that TOR modulates brassinosteroid (BR)-regulated plant growth and stress-response pathways. Active TOR was required for full BR-mediated growth in Arabidopsis thaliana. Autophagy was constitutively up-regulated upon blocking BR biosynthesis or signaling, and down-regulated by increasing the activity of the BR pathway. BIN2 (brassinosteroid-insensitive 2) kinase, a GSK3-like kinase functioning as a negative regulator in BR signaling, directly phosphorylated RAPTOR1B (regulatory-associated protein of TOR 1B), a substrate-recruiting subunit in the TOR complex, at a conserved serine residue within a typical BIN2 phosphorylation motif. Mutation of RAPTOR1B serine 916 to alanine, to block phosphorylation by BIN2, repressed autophagy and increased phosphorylation of the TOR substrate ATG13a (autophagy-related protein 13a). By contrast, this mutation had only a limited effect on growth. We present a model in which RAPTOR1B is phosphorylated and inhibited by BIN2 when BRs are absent, activating the autophagy pathway. When BRs signal and inhibit BIN2, RAPTOR1B is thus less inhibited by BIN2 phosphorylation. This leads to increased TOR activity and ATG13a phosphorylation, and decreased autophagy activity. Our studies define a new mechanism by which coordination between BR and TOR signaling pathways helps to maintain the balance between plant growth and stress responses.

FERONIA functions through Target of Rapamycin (TOR) to negatively regulate autophagy.

FERONIA (FER) receptor kinase plays versatile roles in plant growth and development, biotic and abiotic stress responses, and reproduction. Autophagy is a conserved cellular recycling process that is critical for balancing plant growth and stress responses. Target of Rapamycin (TOR) has been shown to be a master regulator of autophagy. Our previous multi-omics analysis with loss-of-function fer-4 mutant implicated that FER functions in the autophagy pathway. We further demonstrated here that the fer-4 mutant displayed constitutive autophagy, and FER is required for TOR kinase activity measured by S6K1 phosphorylation and by root growth inhibition assay to TOR kinase inhibitor AZD8055. Taken together, our study provides a previously unknown mechanism by which FER functions through TOR to negatively regulate autophagy.

Integration of multi-omics data reveals interplay between brassinosteroid and Target of Rapamycin Complex signaling in Arabidopsis.

Brassinosteroids (BRs) and Target of Rapamycin Complex (TORC) are two major actors coordinating plant growth and stress responses. Brassinosteroids function through a signaling pathway to extensively regulate gene expression and TORC is known to regulate translation and autophagy. Recent studies have revealed connections between these two pathways, but a system-wide view of their interplay is still missing. We quantified the level of 23 975 transcripts, 11 183 proteins, and 27 887 phosphorylation sites in wild-type Arabidopsis thaliana and in mutants with altered levels of either BRASSINOSTEROID INSENSITIVE 2 (BIN2) or REGULATORY ASSOCIATED PROTEIN OF TOR 1B (RAPTOR1B), two key players in BR and TORC signaling, respectively. We found that perturbation of BIN2 or RAPTOR1B levels affects a common set of gene-products involved in growth and stress responses. Furthermore, we used the multi-omic data to reconstruct an integrated signaling network. We screened 41 candidate genes identified from the reconstructed network and found that loss of function mutants of many of these proteins led to an altered BR response and/or modulated autophagy activity. Altogether, these results establish a predictive network that defines different layers of molecular interactions between BR- or TORC-regulated growth and autophagy.

Integrated omics reveal novel functions and underlying mechanisms of the receptor kinase FERONIA in Arabidopsis thaliana.

The receptor kinase FERONIA (FER) is a versatile regulator of plant growth and development, biotic and abiotic stress responses, and reproduction. To gain new insights into the molecular interplay of these processes and to identify new FER functions, we carried out quantitative transcriptome, proteome, and phosphoproteome profiling of Arabidopsis (Arabidopsis thaliana) wild-type and fer-4 loss-of-function mutant plants. Gene ontology terms for phytohormone signaling, abiotic stress, and biotic stress were significantly enriched among differentially expressed transcripts, differentially abundant proteins, and/or misphosphorylated proteins, in agreement with the known roles for FER in these processes. Analysis of multiomics data and subsequent experimental evidence revealed previously unknown functions for FER in endoplasmic reticulum (ER) body formation and glucosinolate biosynthesis. FER functions through the transcription factor NAI1 to mediate ER body formation. FER also negatively regulates indole glucosinolate biosynthesis, partially through NAI1. Furthermore, we found that a group of abscisic acid (ABA)-induced transcription factors is hypophosphorylated in the fer-4 mutant and demonstrated that FER acts through the transcription factor ABA INSENSITIVE5 (ABI5) to negatively regulate the ABA response during cotyledon greening. Our integrated omics study, therefore, reveals novel functions for FER and provides new insights into the underlying mechanisms of FER function.

SGT1-Specific Domain Mutations Impair Interactions with the Barley MLA6 Immune Receptor in Association with Loss of NLR Protein.

The Mla (Mildew resistance locus a) of barley (Hordeum vulgare L.) is an effective model for cereal immunity against fungal pathogens. Like many resistance proteins, variants of the MLA coiled-coil nucleotide-binding leucine-rich repeat (CC-NLR) receptor often require the HRS complex (HSP90, RAR1, and SGT1) to function. However, functional analysis of Sgt1 has been particularly difficult, as deletions are often lethal. Recently, we identified rar3 (required for Mla6 resistance 3), an in-frame Sgt1ΔKL308-309 mutation in the SGT1-specific domain, that alters resistance conferred by MLA but without lethality. Here, we use autoactive MLA6 and recombinant yeast-two-hybrid strains with stably integrated HvRar1 and HvHsp90 to determine that this mutation weakens but does not entirely disrupt the interaction between SGT1 and MLA. This causes a concomitant reduction in MLA6 protein accumulation below the apparent threshold required for effective resistance. The ΔKL308-309 deletion had a lesser effect on intramolecular interactions than alanine or arginine substitutions, and MLA variants that display diminished interactions with SGT1 appear to be disproportionately affected by the SGT1ΔKL308-309 mutation. We hypothesize that those dimeric plant CC-NLRs that appear unaffected by Sgt1 silencing are those with the strongest intermolecular interactions with it. Combining our data with recent work in CC-NLRs, we propose a cyclical model of the MLA-HRS resistosome interactions.[Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2022.

Temporal and spatial auxin responsive networks in maize primary roots.

Auxin is a key regulator of root morphogenesis across angiosperms. To better understand auxin-regulated networks underlying maize root development, we have characterized auxin-responsive transcription across two time points (30 and 120 min) and four regions of the primary root: the meristematic zone, elongation zone, cortex and stele. Hundreds of auxin-regulated genes involved in diverse biological processes were quantified in these different root regions. In general, most auxin-regulated genes are region unique and are predominantly observed in differentiated tissues compared with the root meristem. Auxin gene regulatory networks were reconstructed with these data to identify key transcription factors that may underlie auxin responses in maize roots. Additionally, Auxin-Response Factor subnetworks were generated to identify target genes that exhibit tissue or temporal specificity in response to auxin. These networks describe novel molecular connections underlying maize root development and provide a foundation for functional genomic studies in a key crop.

Integrated omics networks reveal the temporal signaling events of brassinosteroid response in Arabidopsis.

Brassinosteroids (BRs) are plant steroid hormones that regulate cell division and stress response. Here we use a systems biology approach to integrate multi-omic datasets and unravel the molecular signaling events of BR response in Arabidopsis. We profile the levels of 26,669 transcripts, 9,533 protein groups, and 26,617 phosphorylation sites from Arabidopsis seedlings treated with brassinolide (BL) for six different lengths of time. We then construct a network inference pipeline called Spatiotemporal Clustering and Inference of Omics Networks (SC-ION) to integrate these data. We use our network predictions to identify putative phosphorylation sites on BES1 and experimentally validate their importance. Additionally, we identify BRONTOSAURUS (BRON) as a transcription factor that regulates cell division, and we show that BRON expression is modulated by BR-responsive kinases and transcription factors. This work demonstrates the power of integrative network analysis applied to multi-omic data and provides fundamental insights into the molecular signaling events occurring during BR response.

The F-box E3 ubiquitin ligase BAF1 mediates the degradation of the brassinosteroid-activated transcription factor BES1 through selective autophagy in Arabidopsis.

Brassinosteroids (BRs) regulate plant growth, development, and stress responses by activating the core transcription factor BRI1-EMS-SUPPRESSOR1 (BES1), whose degradation occurs through the proteasome and autophagy pathways. The E3 ubiquitin ligase(s) that modify BES1 for autophagy-mediated degradation remain to be fully defined. Here, we identified an F-box family E3 ubiquitin ligase named BES1-ASSOCIATED F-BOX1 (BAF1) in Arabidopsis thaliana. BAF1 interacts with BES1 and mediates its ubiquitination and degradation. Our genetic data demonstrated that BAF1 inhibits BR signaling in a BES1-dependent manner. Moreover, BAF1 targets BES1 for autophagic degradation in a selective manner. BAF1-triggered selective autophagy of BES1 depends on the ubiquitin binding receptor DOMINANT SUPPRESSOR OF KAR2 (DSK2). Sucrose starvation-induced selective autophagy of BES1, but not bulk autophagy, was significantly compromised in baf1 mutant and BAF1-ΔF (BAF1 F-box decoy) overexpression plants, but clearly increased by BAF1 overexpression. The baf1 and BAF1-ΔF overexpression plants had increased BR-regulated growth but were sensitive to long-term sucrose starvation, while BAF1 overexpression plants had decreased BR-regulated growth but were highly tolerant of sucrose starvation. Our results not only established BAF1 as an E3 ubiquitin ligase that targets BES1 for degradation through selective autophagy pathway, but also revealed a mechanism for plants to reduce growth during sucrose starvation.

qTeller: a tool for comparative multi-genomic gene expression analysis.

MOTIVATION: Over the last decade, RNA-Seq whole-genome sequencing has become a widely used method for measuring and understanding transcriptome-level changes in gene expression. Since RNA-Seq is relatively inexpensive, it can be used on multiple genomes to evaluate gene expression across many different conditions, tissues and cell types. Although many tools exist to map and compare RNA-Seq at the genomics level, few web-based tools are dedicated to making data generated for individual genomic analysis accessible and reusable at a gene-level scale for comparative analysis between genes, across different genomes and meta-analyses. RESULTS: To address this challenge, we revamped the comparative gene expression tool qTeller to take advantage of the growing number of public RNA-Seq datasets. qTeller allows users to evaluate gene expression data in a defined genomic interval and also perform two-gene comparisons across multiple user-chosen tissues. Though previously unpublished, qTeller has been cited extensively in the scientific literature, demonstrating its importance to researchers. Our new version of qTeller now supports multiple genomes for intergenomic comparisons, and includes capabilities for both mRNA and protein abundance datasets. Other new features include support for additional data formats, modernized interface and back-end database and an optimized framework for adoption by other organisms' databases. AVAILABILITY AND IMPLEMENTATION: The source code for qTeller is open-source and available through GitHub (https://github.com/Maize-Genetics-and-Genomics-Database/qTeller). A maize instance of qTeller is available at the Maize Genetics and Genomics database (MaizeGDB) (https://qteller.maizegdb.org/), where we have mapped over 200 unique datasets from GenBank across 27 maize genomes. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Robotic Assay for Drought (RoAD): an automated phenotyping system for brassinosteroid and drought responses.

Brassinosteroids (BRs) are a group of plant steroid hormones involved in regulating growth, development, and stress responses. Many components of the BR pathway have previously been identified and characterized. However, BR phenotyping experiments are typically performed in a low-throughput manner, such as on Petri plates. Additionally, the BR pathway affects drought responses, but drought experiments are time consuming and difficult to control. To mitigate these issues and increase throughput, we developed the Robotic Assay for Drought (RoAD) system to perform BR and drought response experiments in soil-grown Arabidopsis plants. RoAD is equipped with a robotic arm, a rover, a bench scale, a precisely controlled watering system, an RGB camera, and a laser profilometer. It performs daily weighing, watering, and imaging tasks and is capable of administering BR response assays by watering plants with Propiconazole (PCZ), a BR biosynthesis inhibitor. We developed image processing algorithms for both plant segmentation and phenotypic trait extraction to accurately measure traits including plant area, plant volume, leaf length, and leaf width. We then applied machine learning algorithms that utilize the extracted phenotypic parameters to identify image-derived traits that can distinguish control, drought-treated, and PCZ-treated plants. We carried out PCZ and drought experiments on a set of BR mutants and Arabidopsis accessions with altered BR responses. Finally, we extended the RoAD assays to perform BR response assays using PCZ in Zea mays (maize) plants. This study establishes an automated and non-invasive robotic imaging system as a tool to accurately measure morphological and growth-related traits of Arabidopsis and maize plants in 3D, providing insights into the BR-mediated control of plant growth and stress responses.

slim shady is a novel allele of PHYTOCHROME B present in the T-DNA line SALK_015201.

Auxin is a hormone that is required for hypocotyl elongation during seedling development. In response to auxin, rapid changes in transcript and protein abundance occur in hypocotyls, and some auxin responsive gene expression is linked to hypocotyl growth. To functionally validate proteomic studies, a reverse genetics screen was performed on mutants in auxin-regulated proteins to identify novel regulators of plant growth. This uncovered a long hypocotyl mutant, which we called slim shady, in an annotated insertion line in IMMUNOREGULATORY RNA-BINDING PROTEIN (IRR). Overexpression of the IRR gene failed to rescue the slim shady phenotype and characterization of a second T-DNA allele of IRR found that it had a wild-type (WT) hypocotyl length. The slim shady mutant has an elevated expression of numerous genes associated with the brassinosteroid-auxin-phytochrome (BAP) regulatory module compared to WT, including transcription factors that regulate brassinosteroid, auxin, and phytochrome pathways. Additionally, slim shady seedlings fail to exhibit a strong transcriptional response to auxin. Using whole genome sequence data and genetic complementation analysis with SALK_015201C, we determined that a novel single nucleotide polymorphism in PHYTOCHROME B was responsible for the slim shady phenotype. This is predicted to induce a frameshift and premature stop codon at leucine 1125, within the histidine kinase-related domain of the carboxy terminus of PHYB, which is required for phytochrome signaling and function. Genetic complementation analyses with phyb-9 confirmed that slim shady is a mutant allele of PHYB. This study advances our understanding of the molecular mechanisms in seedling development, by furthering our understanding of how light signaling is linked to auxin-dependent cell elongation. Furthermore, this study highlights the importance of confirming the genetic identity of research material before attributing phenotypes to known mutations sourced from T-DNA stocks.

Advances in functional proteomics to study plant-pathogen interactions.

Pathogen infection triggers complex signaling networks in plant cells that ultimately result in either susceptibility or resistance. We have made substantial progress in dissecting many of these signaling events, and it is becoming clear that changes in proteome composition and protein activity are major drivers of plant-microbe interactions. Here, we highlight different approaches to analyze the functional proteomes of hosts and pathogens and discuss how they have been used to further our understanding of plant disease. Global proteome profiling can quantify the dynamics of proteins, posttranslational modifications, and biological pathways that contribute to immune-related outcomes. In addition, emerging techniques such as enzyme activity-based profiling, proximity labeling, and kinase-substrate profiling are being used to dissect biochemical events that operate during infection. Finally, we discuss how these functional approaches can be integrated with other profiling data to gain a mechanistic, systems-level view of plant and pathogen signaling.

Post-Synthetic Reduction of Pectin Methylesterification Causes Morphological Abnormalities and Alterations to Stress Response in Arabidopsis thaliana.

Pectin is a critical component of the plant cell wall, supporting wall biomechanics and contributing to cell wall signaling in response to stress. The plant cell carefully regulates pectin methylesterification with endogenous pectin methylesterases (PMEs) and their inhibitors (PMEIs) to promote growth and protect against pathogens. We expressed Aspergillus nidulans pectin methylesterase (AnPME) in Arabidopsis thaliana plants to determine the impacts of methylesterification status on pectin function. Plants expressing AnPME had a roughly 50% reduction in methylester content compared with control plants. AnPME plants displayed a severe dwarf phenotype, including small, bushy rosettes and shorter roots. This phenotype was caused by a reduction in cell elongation. Cell wall composition was altered in AnPME plants, with significantly more arabinose and significantly less galacturonic acid, suggesting that plants actively monitor and compensate for altered pectin content. Cell walls of AnPME plants were more readily degraded by polygalacturonase (PG) alone but were less susceptible to treatment with a mixture of PG and PME. AnPME plants were insensitive to osmotic stress, and their susceptibility to Botrytis cinerea was comparable to wild type plants despite their compromised cell walls. This is likely due to upregulated expression of defense response genes observed in AnPME plants. These results demonstrate the importance of pectin in both normal growth and development, and in response to biotic and abiotic stresses.